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C\) . Abstract 

A fragmenting relativistic string is widely used for modelling particle production 
via quark-gluon strings formed in hadron inelastic interactions of high energies. In 
this note we focus on motion and fragmentation of relativistic string with non-zero 
transverse separation of its ends and study this scenario as a possible mechanism 
bringing anisotropy into the asimuthal angle disribution of produced particles in 
inelastic interactions of hadrons. 
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1 Introduction. 

A two-stage scenario of particle production through the formation and subsequent break-up 
' of quark-gluon strings via Schwinger mechanism is a common picture for inelastic hadronic 
interactions [HE]- A widely used approximation for quark-gluon string is the classical rela- 
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tivistic string with break-ups described by the area law [3|. General picture used in a number 
of models of particle production via quark-gluon string decay (e.g. [U [5]) assumes that a 
quark-gluon string starts stretching from a single point in space though may still have a com- 
plicated structure with kinks in the middle. However, this picture may yet capture not all 
of the features which may be important for some of the observables. Namely, constituents of 
interacting hadrons (valence quarks and diquarks, sea qq pairs) are at different positions in 
the transverse plane during the interaction with a typical separation determined by the value 
of inelastic cross section, that is, around 1 fm, which is considerably less than typical string 
transverse size. 

Another hint comes from the recent studies [Bj, where it was shown that asimuthal 
anisotropy in p t distribution is observed in Gribov-Regge theory for particles produced via 
single Pomeron exchange if the constituents of the incoming hadrons or nuclei, which the 
Pomeron is attached to, have different positions in the transverse plane with the preferred 
in-plane (i.e. parallel to the overall impact parameter) orientation. In hadron-hadron and 
nucleus-nucleus collisions this in-plane orientation appears to be prevailing due to the gradi- 
ents of the optical thicknesses of the colliding systems. At the same time the cut Pomeron 
is usually identified with a couple of quark-gluon strings pQ. So, in present study we address 
a question whether fragmentation of a classical relativistic string with a certain transverse 
separation of its ends could produce asimuthal anisotropy in transverse momentum distri- 
bution of the fragmentation products. We consider motion and fragmentation of a classical 
relativistic string with initial transverse separation of its ends of the order of b = 1 fm. 

The structure of the paper is as follows. In two sections following introduction we review 
generalities of string motion which are well-known and make some general notes which are 
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relevant for our subsequent numerical calculations. The forth section is devoted to motion 
of a relativistic string starting from some two special initial conditions. In the fifth section 
we describe a procedure which we use for string fragmentation in our Monte-Carlo code. 
Numerical results for two types of initial conditions described in section 4 are presented in 
subsequent section. Then follows a conclusion and an outlook. 



2 Relativistic string, equation of motion. 

We start with reviewing the basics of string motion and fragmentation. A detailed derivation 
of the equations of motion starting from the string Lagrangian can be found elsewhere (see 
e.g. [3] or [5] and references therein). Here we just list the results which we shall use in 
subsequent consideration. 

In general, coordinates of the points of the string x M are functions of two parameters, the 
timelike, r £ (—00,00) and the spacelike a (which takes values in a limited interval). To 
constrain freedom in reparameterization one imposes a specific gauge on the x^ as functions 
of a and r. For clarity and transparent physical interpretation of r and a the orthonor- 
mal gauge is convenient, together with so-called lab frame parameterization. The lab frame 
parameterization corresponds to identifying time-like parameter r with time: 

X° = T = t. 

Using this, the orthonormal gauge reads: 

x-x' = 0, m 
v x 2 + x /2 = 1. [L) 

<9x <9x 

where x = — and x' = — — . In this gauge string's equations of motion are just wave equations 
at Oct 

x - x" = (2) 

with boundary conditions 

x'l =0. (3) 

l<T— U,(T max V / 

Energy and momentum of a small piece of a string are 

dp^ = >cx^(t,a)da or dE = xda, dp = xrkda. (4) 

where x stands for string tension. So, parameter a itself has a transparent physical interpre- 
tation, as the conserved energy and momentum of the string are 
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E = x / da = xa max ; (5) 
Jo 

p = x / xdcr, (6) 



JO 

Hence, in the lab frame parameterization a £ (0, E/ x). 

Solution of the equations of motion (EJ) which satisfies boundary condition (jHJ) can be 

expressed as 

*=\[y{t + a)+y{t-a)l (7) 



where y(£), a function of single variable, is obviously the trajectory of one of the string's end 
points. This trajectory is usually called directrix. In terms of the directrix, gauge conditions 
(fTJ are equivalent to 

y' 2 = I- 

Momentum conservation ([6]) implies that directrix is "periodic" : 



This means, we need to know the trajectory of one of the string's ends on a finite time interval 
to be able to fully describe string's motion. At the same time a full picture of the relativistic 
string motion is unambigously defined, if we know initital conditions, that is, momenta and 
arrangement of string at some fixed moment. 

We devote the next section to linking these two points. 

3 Initial conditions and directrix recovery. 

As it has already been mentioned in the introduction, to the best of our knowledge, most 
applicatons use initial conditions for the relativistic string as a source of particles which 
assume that string is stretched from a single point in space. In this case, subsequent string 
motion is fully described once we know momenta of string's ends and kinks at the time of its 
creation. However, it seems natural to consider situations with non-zero length of a string at 
the initial stage of the evolution (non-zero transverse separation for the string's ends) with 
some momentum distribution along the string. One special case of motion of a relativistic 
string with non-zero ends' separation was considered in [8] . Approximation used there assumes 
string's ends moving at constant velocity v < c which physically implies their masses being 
infinite. Contrary, here we shall consider a string with massless ends as a usual approximation 
used in models of hadron production in inelastic interactions. 

We proceed with the technical issue of recovering directrix from coordinates and momenta 
of strings pieces at a given time to in a way which can be used in numerical simulations. Let 
us assume that for this time moment t we are provided information on string arrangement 
in space and momentum distribution along the string length I, x(Z) and These initial 

distributions cannot be arbitrary, the restrictions can be read out from the gauge conditions 
(PQ). One can fulfill gauge conditions having rather x' ^ or x' = 0. The first opportunity 
corresponds to a smooth distribution of -£ with ^ _L -jr. The second opportunity, x' = 
for a certain range of parameter a, implies |x| = 1 with arbitrary direction for x for that a 
range. This corresponds to a 'kink' situation where part of a string of exactly zero length in 
space carries a certain finite momentum leading to a delta-function-like contribution to the 
dp/dl. 

The second of the gauge equations gives a relation expressing increment of the parameter 
a along the infinitezimal piece of a string in terms of its length dl = \dx.\ and momentum dp: 



which upon integration provides dependence /(a). This gives x' and x as functions of cr for a 
given time to: 





da = dE/x = ^Jdp 2 + x 2 dl 2 / 
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for smooth part and kinks respectively. 
At the same time, as follows from (17j) , 



Ato, v) = l [y'(*o + <t)- y'(t - a)} ; (11) 

x(to, <r) = - [y'(* Q + a) + y'(t - a)] . (12) 

So 

y'(*o + <r) =x(t ,cr) + x'(t ,a), , . 

y'(t - ( x)=x(to, ( T)-x'(to,cr). Udj 

Upon integrating ( fl~3l) and taking into account that y(to) — x (to, cr = 0) is a position of one 
of the string's ends, one recovers the directrix on the interval t G (to — 0"max,to + cr max ), and, 
hence, at any point due to its periodicity. These equations can be also written in a form 
convenient for numerical applications, in terms of momenta and lengths of infinitesimal string 
pieces: 

dy(t + a) = y(t + a + da) - y(t +&) = {dp/ x + dx), 
dy(t - cr) = y(t - a -da) - y(t - a) = -(dp/x- dhc). 



4 Examples of motion. 

To illustrate the technique described in the previous section, we reconstruct directrix for two 
simple model initial conditions. These are in some sence extreme cases which implement 
two opposite possibilities, namely a) rapidity gap with large transverse separation between 
partons at strings ends and b) smooth distribution of momentum carried by string pieces with 
linear relation between position of a piece in the transverse plane and its rapidity y. 

We start with the case a). Consider a string with massless ends separated by impact 
parameter b along the 22-axis connected by a string piece being initially at rest. Ends of the 
string carry momenta p and — p along the £3 as shown in fig. [H We call this 'kink-type' initial 
conditions in what follows. In this example total energy of the string is E = 2p + Kb, and 
a G (0, 2p/x + b). The motion of a string with this kind of initial conditions will evidently be 
two-dimensional with trajectory lying in (^2,^3) plane. 

To recover directrix we'll start with drawing a trajectory of the upper endpoint of the 
string for (0, <j max ) time interval, that is using first line of (TH|) - At zero time the endpoint 
carries momentum p with length attributed to endpoint being zero. So for a G (0,p/x) we 
should use (1141) setting c?x = while dp/x differential is non-zero and parallel to X3 axis. 
For a G (p/x,p/x+ b), we have only the second one of the two differentials in the directrix 
increment (fT4j) which goes against x?,- For a G (p/x + b, 2pj k + b) again the first differential 
plays and the increment goes against x%. Directrix pieces for t G (— cr max ,0) are constructed 
in exactly the same manner making use of the second line of (|14j) . 
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Figure 1: 'Kink-type' initial conditions. 



Finally, combining the differentials which is straightforward, we find that the directrix in 
this case represents a rectangular box with argument a changing linearly along the rectangle 
lines. This means that the ends are moving along these lines with the speed of light. Fo lowing 
the prescription (jTJ) we find position of the string at arbitrary moment. The directrix is shown 
in fig. [2]a together with string positions at some selected time moments. Note that a transverse 
part does not disappear instantly and that parts of the strings which are at 45 degrees with 
£3 axis also carry momentum which has a component along x^- So, if fragmentation of string 
starts and two subsequent breaks occur within the inner part cutting off the string piece, it 
will gain an additional contribution to the momentum along x 2 - 
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Figure 2: Examples of string motion with different initial conditions: (a) - 'kink-type', (b) - 
'invariant-type'. For both cases impact parameter b = 0.8 fm, sring cm. energy y/s = 
7.4 GeV, string tension x = 1 GeV/fm. 



For a string piece with longitudinal momentum dp 3 and energy dE one can find its rapidity 
according to usual definition 

1 dE + dps 



2 dE - dp, 
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or in terms of piece velocity: 



» = ilni±^, (15) 
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the latter definition is also suitable for a point on a string except the kinks. 

In the example described above, 'kink' case a), strings ends were connected by a piece of 
a string with exactly zero rapidity in the lab frame. For the second example, case b), we take 
a prescription for momentum distribution along the string which does not lead to appearing 
of some distinguished rapidity value and which is invariant with respect to the boosts along 
the X3 axis. The possible choice which satisfies this condition is to assume that string is still 
a straight line in transverse plane, with rapidity, defined according to ( Tl5j) . changing linearly 
with the transverse position of the point. We denote the proportionality coefficient by a. 
This gives for a momentum and energy of a small string piece of length db: 

dp 3 = ±xdbsmh(a\b — b c \ — y c ). (16) 



This leads to 

dE = xdb cosh(a\h — b c | — y c ) 



(17) 



due to (jUJ). Here b c = (b max + b min )/2 denotes position of the string geometrical center, y c is 
the relative rapidity of the string's center of mass and lab frames while a\h — b c | corresponds 
to the rapidity of the piece with transverse position b in the string center of mass frame. 
Once the string energy in its cm. frame yfs and transverse separation of ends are given, the 
coefficient a can be extracted from 

| b m ax — b m ; n | /2 

\fs = -Ec.m. = J dbcosh(ab) = — sinh ■ max ^ — (18) 

|b max b mm |/2 

e.g. solving it iteratively. 

If the same string is considered in another reference frame boosted along x% with respect 
to the initial one, rapidity attributed to the string point has again a linear dependence on its 
transverse position with the same proportionality coefficient a. The only thing which has to 
be changed in ([16]) is the value of y c . In this sence this prescription for momentum distribution 
along the string is boost-invariant, so in what follows we refer to is as 'invariant-type' initial 
conditions. The distribution of longitudinal momentum (Tl6|) in this example coincides with 
what is found in [5] for the string with infinite masses on its ends at zero time. 

It is also possible to write down an intermediate variant which has both mentioned types 
of initial momentum distributions as limiting cases. Namely one does so setting 

dps = ±A>tdb sinh(a\h — b c |) (19) 

in the string center of mass frame. A is an apriori fixed constant and a is defined to fit a 
particular mass of the string; limits A — > oo and A = 1 give the 'kink' and 'invariant' cases 
respectively. However it is only for A = 1 that the argument of "sinh" function has a meaning 
of rapidity of an infinitesimal string piece and that transformations to other reference frames 
imply just adding rapidity difference to the argument without changing "sinh" law. 

The directrix for the 'invariant' case and string position at any time can be obtained in the 
same fashion as in the previous example. Let us write down the directrix for parameterization 
defined by a and impact parameter in its center of mass frame. Put center of coordinate 
system to the center of the string. Projections of the directrix increments ffl3l) on X2 and x% 
connected with the piece of a string of length dx-i (string lies along X2 at zero time) are 

dy 2 = dx 2 ; (20) 
dy 3 = dp 3 /x = dx 2 sinh(ax 2 ). (21) 

Integration from X2 max downto 22 gives parametric dependence of the directrix on the trans- 
verse coordinate: 

2/2 = x 2 ; (22) 

2/3 = — (coshao^ - coshaa^max), (23) 
a 

where x 2 takes values in x 2 G (— x 2max , x 2max ). The last thing to note is that, the argument 
of the directrix is equal to its length measured from the point at which we start the directrix 
construction due to the condition |y'| = 1 . 

An example of this type of directrix together with string positions at some time moments 
is presented in the figure [2] b. This motion is again two-dimensional and should lead to 
asimuthally anisotropic distribution of fragments on decay. The motion in (22,^3) plane will 
again give additional contribution to particle momenta along the 22 axis. 



As the induced transverse momentum strongly depends on the transverse velocity of the 
string parts, the flow-like anisotropy for both cases considered should be sensitive to both the 
impact parameter and center of mass energy of the initial string, namely growing with impact 
parameter and decreasing with cm. energy. 

5 Fragmentation model. 

To estimate induced anisotropy in asimuthal distribution of the produced particles we use the 
off-shell resonance model (AMOR) [3J with simplifications (the alternative could be effective 
Field- Feynman approach [10]). The AMOR model follows the fragmentation process as it is 
seen from the particular reference frame. At each point of the string a break probability exists. 
A probability of breaking within a string piece within certain time interval is proportional 
to the surface element of the worldsheet swept by the piece in space-time. That is break 
probability within string piece corresponding to da interval (g stands for metric on the string 
worldsheet) is 

— = kdet(-g)da = kda\x'\Vl^i^ (24) 

Cl/L 

with k being some phenomeno logical constant. In this picture (for a reasonable set of param- 
eters k and string tension x giving 1 4- 2 particles per unit rapidity) the string first breaks into 
the so-called primary fragments with mass M 2 ~ 0.5 GeV 2 . As breaking of these fragments 
is an essentially quantum process due their low mass, they are considered as resonances off 
the mass shell and decay (if to neglect spin effects) isotropically in their lab frame producing 
observed mesons (two pions in most cases) [5], [7j . 

The equation ff24l has a transparent physical interpretation, since this means the probabil- 
ity of breaking within some string piece is proportional to its length times the Lorentz-factor 
of the piece. This means in particular, that for two model cases considered in the previous sec- 
tion we shall have different distributions of break points in space and time and consequently 
different asimuthal anisotropy patterns for produced particles. 

Following prescription ( 124|) we make a MC code for string breaking which explicitly ob- 
serves the string dynamics. The basic scheme of the code is as follows. 

1. After directrix is reconstructed numerically from the initial conditions as described in 
previous sections, time evolution of the string arrangement in space is done step by step. On 
each step decay probability is computed by integration of (12"4"I) along the string length, and, 
according to the obtained probability, it is decided whether string will break at this time step 
or not. 

2. If the string breaks at this step, the break point is defined, again making use of (l24"j) . It 
is assumed, that massless quark-antiquark pairs emerge exactly in the break point and have 
zero momentum, so the energy is conserved. This assumption is the main simplification of 
our fragmentation model. 

3. After the breakup point is generated, the two resulting string pieces are considered as 
independent, and are treated in the same fashion as the parent one unless the mass of the 
fragment is below some cutoff m cu t- 

4. If the mass of a newly produced fragment is below the cutoff mass m cut , it decays 
isotropically in its rest frame into two 7r-mesons. This decay is the main source of transverse 
momentum of the final particles in the model (aditional contribution comes from the clas- 
sical string motion), so the transverse momentum distribution is closely linked to the mass 
distribution of the primary fragments. 



Figure 3: Rapidity densities (left) and transverse momentum distributions (right) in our frag- 
mentation scheme with 'invariant-type' initial conditions for different values of impact pa- 
rameter, string cm. energies and cutoff mass m cu t. 



We also do not take into account different flavour and baryonic content of the string decay 
products. Though the model is simplistic, to our mind it keeps the main feature relevant 
for our purpose, which is evaluation of asimuthal anisotropy in momentum distribution of 
fragmentation products. Namely, final particles obtain a random isotropic transverse mo- 
mentum of the order of 7r-meson mass (depends on cutoff m cut ) and an additional 'in-plane' 
contribution to the momentum, which comes from the motion of a string. 

The parameters used in the MC code are string tension x, breaking probability k (see 
(TJl])) and cutoff mass m cut . 

6 Numerical results. 

In our numerical estimates we fix string tension and decay parameter at the level of 



and cutoff mass at m cut = lOm^ as exact fitting of experimental data on multiplicity and 
transverse momentum goes beyond the scope of this work. This set of parameters pro- 
vides reasonable rapidity distribution of average multiplicity with plato at zero rapidity and 
dN/dy\ y=0 ~ 2 (which corresponds to cOVcharged/ dy\ y=0 ~ 1.25) and transverse momentum 
distribution with (p t ) = 0.28 GeV almost independent on the total energy of the initial string, 
see fig. El The transverse momentum distribution drops to zero at approximately 0.5m cut GeV 
as seen from the right panel of fig. [31 which is of course a consequence of the cutoff mass which 
we have in our toy fragmentation model. As one can see, impact parameter b has a minor 
influence on these distributions. 

Asimuthal anisotropy of the transverse momentum distribution is usually described by the 
Fourier coefficients in the decomposition of asimuthal spectrum [TT1 [T2"]: 



where angle is between the direction of transverse momentum and reaction plane which 
contains beam axis and impact parameter. So the flow coefficients v n , provided the direction 
of impact parameter is known, can be evaluated as 



x = 1 GeV/fm, k = 1 GeV/fm 2 




n=l 



(25) 




(26) 
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Figure 4: Flow coefficients v i (upper plots) and v 2 (lower plots) for a decay of a single string 
with 'kink-type' initial conditions. Left: cm. string energy y/s = 40 GeV; right: cm. string 
energy ^/s = 200 GeV. Solid line - impact parameter b = 0.5 fm; dashed line -6 = 1 fm. 
Triangles denote edges of the rapidity distribution defined by ?/ max = In y/s/m n 



with average taken over all particles in the sample. The first V\ and the second v 2 coefficients 
are known as directed and elliptic flow respectively. 

To explore parametric dependence of the flow induced by the two-dimensional string mo- 
tion we perform computations of flow coefficients V\ and v 2 for two types of initial conditions 
described in previous sections for different transverse separation between strings ends to which 
we attribute a role of impact parameter b and different string invariant masses (total energy 
in its cm. frame). We dispose the impact parameter along the x 2 axis with initial momenta 
along £3 and view string fragmentation from its center of mass frame. Hence, to evaluate flow 
coefficients integrated over p t , which we do as a function of rapidity, we compute (cos0) (for 
v\) and (cos 20 (for v 2 ), where <ft is the angle between transverse momentum of a particle and 
x 2 axis which direction coinsides with impact parameter b in our case. 

Results of the computation are presented in two sets of figures. As observed, kink-type and 
invariant-type initial conditions produce different flow patterns. Some qualitative comments 
can be given on the basis of string motion depicted in fig. [2j As suggested in [9], the breakup 
points are causally disconnected, which means that for (not too large) primary fragment 
their transverse momentum closely follows velocity of the adjacent breakpoints which cut this 
fragment away In the 'kink-type' case, fragmentation starts in the region of a string which 
is stretched at zero time between two kinks and which has zero rapidity. The first break 
produces two pieces which have one of their ends slow. Subsequent breakpoints appear on the 
parts of the string which make 45° with the 23 axis cutting off fragments from the slow ends 
of the two secondary strings. This populates rapidity distribution starting from the central 
region. In the beginning fragments with larger lengths are cut off, which hence carry larger 
momentum along x 2 . So anisotropy in this case is stronger at central rapidity. In the case of 
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Figure 5: The same as previous figure for 'invariant-type' conditions. 



'invariant' type initial conditions in-plane transverse velocities of the inner parts of the string 
are much smaller than for the 'kink' case. However ends of the string acquire transverse 
velocities which leads to much more pronounced flow at forward and backward rapidities. 
This is illustrated in figure [6] showing the in-plane transverse velocity distribution of primary 
fragments (produced in the model via classical string fragmentation as described in previous 
section) prior to final decay into two particles. 
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Figure 6: Rapidity and in-plane transverse velocity /3± = v ™& men j c distribution of primary 
fragments: 2D plot, showing number of counts as a function of transverse velocity and rapidity 
N(/3 t ,y). Surface of the box on the plot is proportionnal to bin content. Left: 'kink-type' 
initial conditions; right: 'invariant-type' initial conditions, string cm. energy i/i = 40 GeV. 
Both cases are for impact parameter 6=1 fm. 

It is remarkable that the sign of the directed flow coefficient v\ in both cases is opposite 
to the 'spectator flow', that is, sign of the directed flow is the same as in the experimentally 
observed directed flow for pi-mesons in Au-Au collisions [13J. 



7 Discussion and concluding remarks. 



String motion and fragmentation started from two initial conditions depicted here represent 
two in some sence extreme cases. For the 'kink-type' initial conditions almost all of the string's 
energy is initially carried by its ends while for the 'invariant-type' the energy is distributed 
in the interior of the string. Hence the value for the elliptic flow at central rapidity obtained 
for the 'kink' case (fig. H]) should be rather viewed upon as the upper estimate for the part of 
the flow originating from the anisotropy of the quark-gluon string decay itself. In general it is 
remarkable that in both scenarios we obtain positive elliptic flow and the sign of the directed 
flow coincides with that experimentally observed. 

What is also remarkable, it is the sensitivity of the rapidity dependence of the flow in the 
model to the initial momentum distribution carried by the interior of the string, as illustrated 
by different rapidity dependencies of the flow depicted in figs. HlEl We note, that this 
statement should be valid for any fragmentation scheme, e.g. Field- Feynman [10], since, as 
mentioned, break points on a string are causally disconnected. 

Application of the developed technique to the interaction of hadrons or nuclei is still not 
straightforward as requires careful treatment of several important features, such as, e.g., fluc- 
tuations in number and invariant masses of strings, produced at the first stage of interaction. 
A careful treatment of initial conditions with paying special attention to longitudinal momen- 
tum distribution in the inner part of a string and its correlation with the transverse position 
is probably the most important part. The fragmentation model itselt may also require future 
refinements, such as implementation of transverse momentum of q-q pair produced at breakup 
of a string, which is for simplicity set zero in this work but is necessary for more realistic p t 
distributions of final particles and accurate study of the p t dependence of the flow. 

These options however lie beyond the scope of the article which has a qualitative character. 
Its aim is to suggest a new mechanism which can give contribution to the flow in the models 
based on quark-gluon string picture of particle production. It also could be important for 
setting up initial conditions for hydrodynamic evolution. So we point out that transverse 
separation of string ends and correlation between rapidity and transverse position of string 
parts can play an important role. 
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